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The Klein-Gordon equation is interpreted in the de Broghe-Bohm manner as a single-particle 
relativistic quantum mechanical equation that defines unique time-like particle trajectories. The 
particle trajectories are determined by the conserved flow of the intrinsic energy density which can 
be derived from the specification of the Klein-Gordon energy-momentum tensor in an Einstein- 
Riemann space. The approach is illustrated by application to the simple single-particle phenomena 
associated with square potentials. 



I. INTRODUCTION 



In a recent paper ||] we outlined a causal trajectory interpretation (in the de Broglie-Bohm sense) of the Klein- 
' Gordon equation. We indicated how to proceed in order to retain the notion of particle trajectories, both for the 
CN • single- and many-particle cases. The basic idea is to extend to the quantum context some results established in the 
context of general relativity for conserved flows of energy-momentum in classical scalar fields. In the classical context 
Edelen [^has shown that there is a natural definition of rest-energy flow and an associated conserved density which 
I are determined respectively by the time-like eigenvectors and eigenvalues of the energy-momentum tensor of the scalar 
7—1 . field itself in its local rest frame. In the context of a single-particle Klein-Gordon equation these time-like eigenvectors 
^-H ' allow the definition of well-behaved particle trajectories. This approach overcomes some of the previous objections to 
. the interpretation of relativistic scalar fields in terms of a particle ontology. 

' In this paper we give further details of our approach and provide an illustration by showing how the simple quantum 
. phenomena associated with single-particle scattering from static square potentials can consistently be interpreted in 
I terms of well-defined individual relativistic boson trajectories. The single-particle picture is consistent and provides 
'■^ ■ no requirement to have recourse to second quantization. (Particle creation can only be discussed in a many-particle 
approach.) We also discuss the Lorentz transformation properties of the trajectories. 



^ ■ II. PARTICLE INTERPRETATIONS OF THE KLEIN-GORDON EQUATION 



de Broglie first proposed a particle interpretation of the Klein-Gordon equation in the period 1926-1927 Q. de 
Broglie proceeded using the polar form of the scalar field, 



■ (/. = i?exp(i5) (1) 

to decompose the wave equation 

Ucj) = -ml(j> (2) 



into a continuity equation 



and a Hamilton-Jacobi equation 



df'iR^d^S) = (3) 



d.Sd^^S -^^ml (4) 
K 

This approach gives a positive definite Lorentz scalar probability density, , but the problem with the use of d^S 
to define the flow lines is that d^S is not always a time- like four vector. Gonsequently the three- velocity, defined (in 
the absence of electromagnetic fields) by 



1 



v^-^ (5) 

dt 

is not necessarily subluminal. This fact was, of course, known to de Broglie |^ but he did not consider it a serious 
pathology. After all, no experimental consequences follow from the implied superluminal motion of the particles and 
the measurable predictions of the quantum theory are recovered and these do not contravene relativistic covariance. 

Another approach based on relativistic hydrodynamics has been developed ( Q], 0) but in this work also 
there is an implied assumption that df^S is a time-like four- vector. Because of the difficulty with space-like flows the 
particle interpretation of the Klein Gordon equation has not generally been accepted. Instead an interpretation based 
on well-defined, extended fields has been developed g, g, (lO|, 

It turns out however that at least for massive bosons there is a possibility to define time- like flows and hence particle 
trajectories that follow the flows and so one objection against a particle formulation of boson fields may be overcome. 
Our aim here is to explore the consequences of taking the idea of boson trajectories seriously. 

Dirac JT^ has shown that in relativistic hydrodynamics a Hamiltonian description can be given with reference to 
an arbitrary space-like hypersurface. The Hamiltonian density H for any field splits into two pieces, each piece being 
proportional to a lapse or shift vector (TV, Ni). The two pieces have a simple physical significance as shown by Schutz 
[Q, namely 

N^^^^^jf^ = {-9)--n^T^,,n''N = eN (6) 

iV.ffl^^ = i~9)h'''n^T,,N, ^ P^N, (7) 

where = —N{^g^'^) is a unit normal to a space-like hypersurface, i^g^") is the four dimensional metric and g*-' is 
the three-dimensional metric, e and Pi are, respectively, the co-ordinate densities of energy and momentum measured 
by an observer at rest on the hypersurface. 

One can define the flow lines for relativistic bosonic fields through the energy- momentum tensor Tj^. In order to 
obtain a Lorentz invariant description one must form true 4- vectors to define the flow of energy momentum p2[ . If 
one defines an arbitrary time-like vector if then Tjfrj'^ is a true four-vector defining the flow of rest energy.[J. This 
flow, however, is dependent on the arbitrary choice of 7]'^ and is not therefore suitable, as it stands, if one wishes to 
ascribe beables to the field itself. 

Edelen has suggested that an intrinsic natural four-vector is provided by the matter field itself through the 
eigenvalue equation 

Tl'W = \W^' (8) 

Edelen shows that, given such a unique W^, an intrinsic rest energy density {p) exists and does not have to be 
introduced ad-hoc into the Einstein theory. The intrinsic energy density satisfies 

{pW'^),^ = (9) 

where ();^ indicates the covariant derivative. Edelen establishes, given the unique W^, that 

T^'" ^ XW^'W + (J^"' (10) 

where a^'^ is any symmetric tensor field with space-like support. Using also (T^'');^ ~ one gets 

{XW),, - rr^-^e^, (11) 

where e^^ is the Born rate of strain tensor. One now looks for a path density S such that its Lie derivative along W^^ 
equals {—g)^a'^'^ef_i^. It follows straightforwardly that 

5= (exp[-^ edr]) ^ (exp[^ edrji-gf^a^'^e^^dr^ (12) 



'^A justification for this approach can be found in M and uA 
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all integrations being performed along the and e = e^^g^'^ . 
In the presence of electromagnetic fields we have 

where is the charge current, F^" is the electromagnetic field tensor and d^^F^'^ is the Lorentz 4-force. The modified 
conservation relation now reads 

{XWn,. = ^^"e^,. + ej^F^^'^W, (14) 

The path density is now defined to be 

S{so) = (exph / " erfr]) ( / (exp[ / edr](-g)^ {(J^" e^, + j^F^^'^W,) dr] 

S will, in general, be path dependent and if the quantities in the integrand are continuous with no singularities can 
be set to zero on an initial space-like hypersurface. Therefore the initial value of p is equal to A and A is a constant 
of the motion along the flow lines. 



One now has a unique scalar density 



or, if one prefers, a unique p given by 



$ = (-g)5A-S' (15) 



$ = (-.g) V (16) 



One can say that \{-g)2W^' is a flux of rest energy density created in amounts equal to a^'^e^^i, per unit geometrical 
volume, cr^'' are the generalised stresses giving rise to the energy flow. S then represents a path-dependent flux which 
combined with the flux of rest energy density gives a conserved flux . can thus be thought of as an intrinsic 

energy density flux. 

Using this approach for the Klein-Gordon field circumvents a number of objections to the very notion of boson 
trajectories, either using the de Broglie approach or using the components of the energy-momentum tensor to define 
flow lines and densities 

Other objections have been raised to single-particle formulations of relativistic quantum mechanics. For example 
there are a number of ways of setting up a Hilbert space for the Klein-Gordon equation ^ but the scalar product 
utilised does not lead to a positive definite probability density. The position operator, also, does not have the usual 
form and it is not clear what measurement procedure would correspond to the measurement of position in the theory 
since the eigenfunctions arc not delta functions and have an infinite range which suggests a violation of causality. 

One way of proceeding has been to use the current which obeys a conservation relation 

d^fd^x^ I fn^d^x (17) 



V 

where V is a four volume and S is an enclosing hypersurface with unit normal n^ . Since is not time-like everywhere, 
even in the case of a restriction to positive energy solutions, j'^n^ does not provide a "probability" density analogous 
to the non-relativistic expression. In the case of trajectory theories one could use a scalar density j'^n^ to give a 
density of crossings of the flow lines, given by j^, across a hypersurface. The usual choice is to use equi-time surfaces 
giving a density of j^rto (with ?io = 1). A Lorentz transformation to a different frame of reference and equi-time 
hypersurface then leads to a density p', given by 

p' = J>"' = P (18) 

One has however privileged one equi-time hypersurface (one does not choose jo in each frame). This deflntion of p 
does not, therefore, depend solely on the held but also, through rt^ on the choice of the space-like hypersurface. (Diirr 



^See for example Schweber Wald An example of a different approach is that of Gitman et al. 
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et al. propose to treat the hypersurface as an additional dynamical variable with its own equation of motion.) 
One should note that the usual Klein-Gordon scalar product, defined via the conservation relation, is independent 
of the choice of surface S and if one considers an infinitesimal tube along the direction of j'^ then putting = jj^ 
one will have |j^|(icr constant along the tube, where da is the cross-sectional area normal to j^. The fact that is 
not everywhere time-like, and hence the surfaces a not always space-like, means that jj^j is not interpretable as a 
probability density of finding a particle in a given position. The density of crossings of flow-lines is however given 
by \jfi\ which is a scalar quantity. In one frame of reference, at one point in space-time, there will be a rest frame 
in which p — in the case of a time-like flow. Although this is not of use in the Klein-Gordon case it suggests 
that, in the Dirac case for example, one might choose the density of crossings as Ij^^l which will be independent of the 
hypersurface and Lorentz invariant. 

Although there may well be a preferential slicing of space-time [|o| and equations of motion associated with that 
slicing, in view of the lack of a theory of such a preference, it seems better to suppose that the choice of space-like 
hypersurface for the preparation of a system in a given state and its subsequent measurement can be left open to the 
experimentalist (in thought, at least). 

We therefore have proposed to take the eigenvalue A of the energy- momentum tensor Tjf, corresponding to the unit 
time-like vector W^^ as the positive definite scalar density. In the Klein-Gordon case we have previously given the 
prescription for finding which always exists and solely depends on the structure of Tjf. X being given by :- 



ml 101" -f 19^09^01 (19) 



A is the rest energy at a point, in the frame of reference provided by [W^] = [1,0,0,0] i.e. the locally stationary 
frame, and X^y—gW'^'■ is a flux of rest energy density p^/—g such that 

ipW");^ = (20) 

It is important to note that the flow lines given by are not geodesies since in the Klein-Gordon case in equation 
(|l^) the (7^'' is non-zero i.e. the Klein-Gordon field is not that of an incoherent fiuid - there are additional state- 
dependent pressure terms of a purely quantum mechanical origin. One sets p equal to A on an arbitrary initial 
space-like hypersurface; P\/—g has been called the intrinsic energy density The intrinsic energy density per unit 
geometrical volume is constant along the fiow lines defined by ^ i.e. \^—gda is constant along the flow lines and 
changes in da are compensated by the ^—g factor, where 



■' =exp{ + I eds ] (21) 



VF9% 

and e = , being the Born rate of strain tensor. 

A is, therefore, a constant of the motion (2). The advantage of using A as a density of crossings is that it is defined 
solely by the field and is a constant along the time- like fiows; one may also consider arbitrary space- like hypersurfaces 
for preparation and measurement. 

One must however note that A is a non-linear observable of the type proposed by Weinberg |Q and Leifer ]2^ ] 
who give different accounts of a non-linear modification of quantum mechanics. Weinberg proposes a scheme that 
would enable one to measure such a non-linear observable. A is clearly not a conventional probability density. In 
justification of such a procedure one notes that in trajectory theories, such as those proposed by de Broglie and Bohm, 
there is no need to use the conventional apparatus of Hermitean operators and collapse to define the theory. One 
has a well-defined physical picture in which the particle follows a trajectory and once a measurement interaction is 
specified the outcome can be calculated for various initial particle co-ordinates. The uncontrollable nature of the 
initial co-ordinates in any given case limits the prediction to statistical statements only. One may therefore propose 
that the actual density of crossing of flow lines corresponds to that of the particles with no necessity to consider 
conventional measurement processes and no need to define a probability density yielding the statistical results if one 
were to make a measurement. One takes a consistent ontological position and accepts the contextualisation of results 
upon measurement. 

The eigenvectors of Tj^ are defined by . In the case of massive spin and spin 1 fields unique time- like eigenvectors 
exist at each point of the Riemannian curved space-time. Transforming to Riemann normal coordinates gives in 
the usual fiat space-time form (minimal coupling with no curvature contributions p6[ . We discuss here only the scalar 
case. 

Writing a solution of the Klein-Gordon equation as 
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(j) = exp[P + iS] 
and setting d^P — P^j, and d^S = one gets 

Ti; = \<j>\^[ml ~ {P^P^ + S"Sc)K + 2|0n(P''P. + 5^5,)] 
The second term can be written as 2|(/)p times 



(22) 



(23) 



(24) 



the eigenvectors of TIf must be a Hnear combination of [P^] and[S'^] (or orthogonal to the term above). Some algebra 
shows that the two solutions are of the form (unnormalised) : 



5^ + exp'' P'^ 



where 



nnhO — 



2Pf^S„ 



(25) 
(26) 



(27) 



These two vectors are easily shown to be orthogonal four vectors using the expression for 
sinhO. One can conclude that one is time-like and the other space- like (or both null). If both eigenvectors were 
space-like then, in some frame of reference, their time components would both vanish giving — —e^^ ~ "fb ^'^^'-'^ 
is not possible with 6 real as given by (27). 

In deriving one has assumed the usual metric with one time-like axis. As a result two further space-like vectors 
and 6^ can be derived such that 



e'P 



and similarily for 6^, since the vectors 5^ 
overall eigenvalue A is given by 



e P^ and S*^ — e 



= 
= 

'Pp define a space-like two dimensional plane. The 



A = yml ± ([P^P^ - S'^S^f + 4[P^5^]2) 

The three- velocity is given (in the absence of electromagnetic fields) by one of the two expressions 

_ _VS^±_e±^VP_ 



(28) 



(29) 



In those regions where the de Broglie three- velocity is less than 1 then 9^5* is a time-like four vector in all frames 
of reference. A covariant condition can then be given for agreement between VdB and v^- Transforming to the rest 
frame at a point and rotating the space axis one gets 

ept _ _ ql p-Opi _ qt 



The second case can be dealt with by setting 6 equal to —9 and P^ equal to — P^ in the first case, so one need only 
consider the first case in detail. If one transforms back to a frame moving with three- velocity equal to tanha (at 
some point)onc finds: 
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= -Pq (e^ + e"^) sinhacosha - Pi (cosh^ ae^^ + sinh^ e^) 
So = Po (cosh^ ae*^ + sinh^ e"") + Pi [e'^ + e^'^) sinhacosha 

If the de Broglie three- velocity is less than one (at the considered point) one sees that 6 must be sufficiently large and 
negative to get 



_Si_ _ 

—So 

The two ways of calculating the velocity will agree where 

Si 



-So 



tanh ( 



< 1 



with 9 sufficiently large and negative (or positive for the other case). In one dimension it can be shown that the + 
signs give the time-like flows. The flows of the field and hence the individual particle trajectories are described by 
a set of world-lines in space-time. Evidently they transform according to the appropriate Lorentz transformation to 
any other inertial frame. The world-lines of the particles are Lorentz invariant in the sense that the same set of events 
will be connected by the transformed world-line in all inertial frames. 



III. SIMPLE EXAMPLE: THE SQUARE POTENTIAL BARRIER 

In order to illustrate our physical model we compare the velocities and trajectories defined using the non-relativistic 
de Broglie-Bohm formulation, de Broglie's approach to the Klein Gordon equation and our definition in the context of 
single-particle scattering from a square potential barrier. The counter-intuitive behaviour which arises in this context 
is often cited as evidence of the inadequacy of single-particle interpretations of relativistic wave equations and of the 
need to proceed to second quantization, but on closer examination it can be seen that no paradoxical behaviour arises 
provided one determines the boundary conditions appropriately. 

Relativistically the potential may be either scalar or electrostatic in origin and we consider both cases here ^. Using 
the relativistic relation 

{hujf = {nkf + ml (30) 
(we take h = c = I). In the regions of zero potential the momentum is given by 

ki = ±\/w2 - (31) 
, The momentum in the potential region depends on the type of potential. For the scalar case we have 

ks = ±^jLU^-{mo + Vf (32) 

whereas for the electrostatic case we have 

ks - ±^ (lu ~ Vf - mo^ (33) 

When dealing with an electrostatic potential step of semi-infinite extent one needs to take the positive square root 
in defining the momentum in the potential region in order to ensure that we have the correct physical boundary 
conditions. This choice avoids the appearance of waves travelling in from positive infinity when E < V + m and hence 
the Klein paradox does not arise. R We also take positive values for w, there is an exactly symmetrical set of solutions 



^The scalar potential was originally introduced to give a containment model for quarks (the so-called bag model) after it was 
found that an electrostatic delta-function barrier could not give containment. 

''The original Dirac theory version of the paradox was effectively resolved in this manner |^ , see also the discussion in ||2^] 
p264. Pair production by strong fields is a reality, but evidently it can not be treated in a single-particle theory restricted to 
the single-particle sector of Fock space. Examples involving pair production can only be treated in many-particle theories. 
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for negative lo but the two sets are disjoint. We shall consider in detail only the barrier case here. The electrostatic 
step behaviour cannot be obtained by extending the barrier width to infinity as one always has waves propagating in 
each direction when E < V + m. 

We define the potential barrier to occupy the region < x < a and to have magnitude F in a frame of reference 
that we label S. In front of the barrier, for a; < the wave function is 

V'(x) = e''" ^ + Re-''" ^ (34) 

whereas in the region of the barrier we have 

^{x) =Ge''"''' +H (35) 

beyond the barrier 

^-(2;) = Je''" ^ (36) 

Solving the boundary conditions at the barrier we find the standard expressions for the reflection and transmission 
ratios. 



\Tf = ^^^^ (37) 



, ,n (k? — kn) siv?k2a , , 

\R\ = T 38 

[kl + klf -{kl-klf cos^{k2a) 

The different cases studied here are distinguished by the prevailing relationship between k and oj. In the scalar 
case k2 becomes imaginary when V > u) — m (and the so-called Klein paradox does not arise for scalar potential 
steps), whereas in the electrostatic case k2 is imaginary only in the interval uj — m <V < lo ^ for which there is 
exponential decay in the barrier. As the barrier height is increased beyond w + to, k2 becomes real and one has a 
transmitted and a reflected wave in the barrier region and there is no confinement. 

In the usual approach the physical meaning of the component wave functions is derived from their associated 
currents, one speaks in terms of an incoming, transmitted and reflected current. But the wave function in front of the 
barrier is, of course, a pure superposition state and so strictly there is just the overall current associated with this pure 
state. Where a plane-wave description is used, the incoming and reflected wave functions overlap over the whole of 
space (x < 0) and so, in the usual approach, for the purpose of interpretation the wave function is tacitly considered 
to be a mixture; only then can the physical situation be described in terms of separate incident and reflected currents. 
The reflected current is given by — fcj \R\^ and the transmitted current by 3?(fc2 |7^|^)- In the pure state the current 
is always positive or zero. As the reflection coefficient increases the current decreases, vanishing as the refiection 
coefficient approaches unity. 

For the scalar potential the variation of the transmission and reflection coefficients for a given energy and over a 
range of barrier heights and widths is similar to that familiar from the non-relativistic case, ^e set m = 1 for the 
calculations presented in this section. As shown in Fig. 1. the transmission decays to zero as the barrier height 
increases to V = lu — m and the wave is completely refiected. For an electrostatic barrier the situation is somewhat 
counter intuitive. For an incident momentum of 0.95 (in our arbitrary units), Fig. 2. shows the behaviour of the 
transmission coefficient (for the electrostatic case) as the magnitude and width of the barrier are varied. At first, as 
expected, the transmission decreases as the barrier potential height or width increases and this continues until the 
region for which V = uj ± m throughout which there is an exponential decay of the wavefunction in the barrier. As 
the magnitude of the barrier potential increases still further, transmission once more increases right up to = 2uj 
where full transmission occurs. Thereafter transmission oscillates with increasing barrier height and width (where in 
the scalar or non-relativistic case the transmission is zero). 

In the non relativistic case the continuity equation for the conserved density p = \tp\'^ is 

|f-V.(p«-) (39) 



5W 
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and this suggests that the velocity be defined by 



P '(Pix,t) 

which for the region in front of the barrier yields 



vs - == M&P^ (40) 



~ 2R{cos{2ki x)) + l + R^ ^^^^ 



At the minima of the denominator the velocity is given by 

ki (1 + R) 

Evidently for the non-relativistic case the velocity has no upper bound and where the density is very small the velocity 
will be very large. The density will develop nodes as i? ^ 1. This case was solved exactly by Takabyasi [|3|. Typical 
trajectories for the incident wave packet case can also be seen in p7|] . 

De Broglie's expression for the velocity in the Klein Gordon case is given in equation (||). For the particular case 
of the barrier this yields 

VdB = 7 — ^—TT^vs (43) 
{uj - V) 

with ki in the expression for vs now given by (p3|), for the region in front of the step. In the limit for which fci is small 
the de Broglie velocity approaches vs- Again there is no upper bound on the velocities defined in this way. As the 
reflection coefficient approaches unity the regions around the minima for which the velocity is superluminal shrink, 
but the magnitude of the velocity in these regions increases rapidly. The regions shrink to a point as the reflection 
coefficient approaches unity. 

The situation is markedly different if one uses the time-like eigenvector of the energy-momentum tensor to define 



velocities and its eigen- value to define the density as described above in equations (29) and (^8|). In this particular 



one-dimensional case, for the definite energy eigenfunctions, the expression for the velocity can be written 

VS + e^VP 

Ve = 44 

uj — V 

This velocity is a weighted combination of the de Broglie velocity and what has been referred to elsewhere as the 
osmotic velocity. The combination always has a magnitude less than 1. 

For this simple case where the wave function in front of the barrier is a superposition of two counter-propagating 
waves with complex weights, it can be shown that the eigenvalue of the energy- momentum tensor A and are in 
phase and related at the maxima and minima by 

A = 2fc2|i?| + (l±|i?|)2 (45) 

In this case then, even for complete reflection when i? = 1, although \4>\^ will go to zero at the minima, A remains 
finite everywhere. f\ If one wants to calculate p at later times then he piccewise constant potential introduces some 
complications since the delta-function forces at the beginning and end of the barrier introduce discontinuities in the 
path density S and one will have to integrate along the paths. We have 

■"'^"'"--^■'•(^-^)"'" ,46, 

For the one-dimensional, time-independent electrostatic barrier this becomes 

j^F.^W'^ = e ifW' - fW') ^ (47) 



^This is similar to the Dirac case where the density can not have nodes 
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which has a singularity at the start and end of the barrier, hence there will be a step discontinuity in S at the beginning 
and end of the barrier. 

For the purposes of illustration we calculate the various densities, trajectories and velocities defined above, in the 
lab frame in which the barrier is stationary, for a variety of cases. 

For a momentum of ki =0.1 (in our arbitrary units in the S frame) the particle has a velocity of one tenth of the 
speed of the light. If the reflection coefficient is also low then the three expressions for the velocity are approximately 
equal, as are the densities and the motion is time-like in all cases. The situation is very different if we consider a large 
reflection coefficient, say R = 0.99. With the relativistic energy given hy E = vLOl the Schrodinger velocity and 
the de Broglie velocity are similar, the relativistic scaling is just vLOl. Both expressions give superluminal motion 
in the interference minima. The eigen-vector velocity remains sub-luminal showing only a small oscillation. As the 
reflection coefficient approaches unity the de Broglie velocity increases without limit at the minima of \(j>\^, whereas 
the energy-momentum velocity approaches zero as A approaches a minima. A graph comparing the de Broglie and 
the energy-momentum velocities and densities is given in Fig. 3. It is clear that even in this " non- relativistic" limit 
the relativistic corrections to the velocity are important in maintaining a sub-luminal velocity in the interference 
minima. That the low-energy limit of the relativistic de Broglie-Bohm velocity and the non-relativistic de Broglie- 
Bohm velocity differ considerably in interference situations was first noted in |Q.The de Broglie trajectories associated 
with this situation are given in Fig. 4., and the eigen-vector trajectories in Fig. 5. The relativistic trajectories at 
low energy differ considerably from those calculated using the non-relativistic (or even the de Broglie) formulation. 
Several authors have suggested using non-relativistic Bohm trajectories to calculate low-energy barrier-tunnelling 
times, however we see that the calculation of tunncUing-times based on the low-energy relativistic trajectories will in 
general give rather different predictions. 

At the relativistic momentum k = 0.95 the differences are more marked for all values of the reflection coefficient. 
Fig. 6 shows a comparison of the velocities and densities defined in the de Broglie theory and in ours for the case of a 
reflection coefficient of 0.7 when E < V ~ m (the normal region). Fig. 7. shows the associated de Broglie trajectories 
and Fig. 8. the associated energy-momentum trajectories. Note that the velocities differ appreciably in the region 
of the interference minima. When V — mKEKV + m the reflection coefficient is unity and the velocities are 
zero. When V > E + m we enter the Klein paradox region. Fig. 9. shows the de Broglie and eigen-vector velocities 
and densities for a reflection coefficient of 0.7 in the anomolous transmission region. Fig. 10. shows the de Broglie 
trajectories associated with the velocities of Fig. 9. Note that the motion is always in the positive direction. Fig. 11. 
shows the eigen-vector trajectories associated with the velocities of Fig. 9. 

IV. LORENTZ COVARIANCE 

In the theory we have given here, based on the eigenvalues and eigenvectors of the energy momentum tensor in 
the local rest frame of the flow, the density is a Lorentz scalar and the flows of energy-momentum are defined by 
four-vectors and are always time-like. The Klein-Gordon wave function, which determines the density and the flow 
evolves unitarily in a covariant manner and so the problems associated with wave-packet collapse in a relativistic 
theory do not arise in the same way (there is no wave packet collapse). If the particle position (or the entire world 
line) is taken as the beable or the element of reality then these are clearly Lorentz covariant. For a given space-time 
wave function $ (x) — <&' {x'), there is a unique world line through a given space-time event. The set of space-time 
points linked by a world line is invariant, but their co-ordinate description in other frames is, of course, different. So 
the situation here is no different than that for a classical particle world-line. We labour this point as doubts have 
been raised concerning the Lorentz covariance of trajectories as deflned in relativistic de Broglie-Bohm theories. 

To expand on this point further, consider a single-particle system prepared in the laboratory or S frame on a 
space-like hypersurface (given by the actual experimental details), typically this will be an equal-time hypersurface, 
although it need not be so. It is important to remember here that all preparations must be referred to some extended 
space-like hypersurface and not to a single point or to a time-like hypersurface. From this initial speciflcation on 
the arbitrary space-like hypersurface the evolution of the wave function can be calculated and the space-time energy- 
momentum flow lines are then determined using the prescription given above. The calculation of the wave function 
on any other arbitrary space-like hypersurface, different from the initial surface, would have to be carried out using a 
type of Tomanaga-Schwinger approach, evolving the initial wave function to the chosen hypersurface. The flow lines 
remain the same, irrespective of the choice of hypersurface, although their co-ordinatisation using alternative frames 
of reference will be different. There is an important distinction to be made between the arbitrary choice of space- like 
surface and an inertial given frame. 
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Consider now the same physical situation as it is described in another inertial frame, E', in motion with respect 
to S. In this frame the system is not described as having been prepared on an equal-time hypersurface, the initial 
hypersurface in E' is just the Lorentz transformed hypersurface defined in E. The relationship between the flow-lines 
in E and E' is simply that of the usual passive re-coordinatisation of world lines which follows on the application of 
the appropriate Lorentz transformation. There is no preferred frame in this approach, although there is, as a matter 
of contingent fact, a particular space-like hypersurface on which the system is actually prepared. 

It is perhaps not surprising that Lorentz covariance follows naturally in the single-particle case where non-local cor- 
relations cannot arise. We consider the many-particle case in a further paper. In the approach that we have described 
here measurements must be treated as particular physical interactions described by introducing the appropriate inter- 
action terms in the Hamiltonian and the appropriate additional apparatus variables. There is no wave packet collapse 
and the wave equation will still yield a covariant process evolving unitarily into different channels each associated 
with a different outcome of the measurement In a further paper we will consider the Klein-Gordon equation for 
a fixed number of particles greater than one and its interpretation in terms of particle trajectories. This will allow 
consideration of Lorentz covariance and nonlocality in a relativistic boson particle trajectory de Broglie-Bohm theory. 



[1] C. Dewdney and G. Horton (1996) in Bohmian Mechanics and Quantum Theory: An appraisal, edited by J.T. Gushing, 
A. Fine and S. Goldstein. Boston Studies in the Philosophy of Science 184. 
D. G. B. Edelen, Nuovo Gimento XXX, 292 (1963). 
L. de Broglie, C.R. Acad. Sci.183, 447 (1926); 185, 380 (1927). 
L. de Broglie, C.R. Acad. Sci. 185, 1118 (1927). 

L. de Broglie, Nonlinear Wave Mechanics (Elsevier, Amsterdam, 1960). 
F. Halbwachs, Th'eorie relativiste des fluides a spin.{ Paris 1960). 
J. P. Vigier and A. Garuccio, Lett. Nuovo Gimento, 30, 57 (1981). 
D. Bohm, Phys. Rev. 85, 166 and 180. 

D. Bohm, B.J. Hiley and P. Kaloyerou, Physics Reports, 144, 6, 349 (1987). 
M. M. Lam and C. Dewdney, Found. Phys. 24, 3 (1994). 

D. Bohm and B. J. Hiley, The Undivided Universe (Routledge, London 1994). 
P. R. Holland, r/ie quantum theory of motion (GUP, Cambridge 1993). 
P.A.M. Dirac, Proc. Roy. Soc. London. A246, 326(1958); A246, 333 (1958). 
S. Schutz, Phys. Rev. D 4, 3559 (1971). 

M. S. Masden, Classical and Quantum Gravity 5, 627 (1988). 

S.S. Schweber, An introduction to Relativistic Quantum Field Theory (Harper and Row, New York 1961). 
R. M. Wald, Quantum theory in curved space-time and black hole thermodynamics (Univ. Chicago Press, Chicago (1994). 
D. M. Gitman and I.V. Tyutin, Classi cal and Quantum Gravity 7, 2131 (1990). 
D. Diirr, S. Goldstein and N. Zanghi, |quant-ph/9801070| ,1999 (unpublished). 
A. Valentini, Ph.D. thesis, London University, (1992). 
H. G. Dosch, J.H.D. Jensen and V.L. Mueller, Phys. Norv. 5, 151 (1971). 
W. Greiner, Relativistic Quantum Mechanics. (Springer- Verlag, Berlin Heidelberg 1990). 
T. Takabayasi, Prog. Theor. Phys. 9,187 (1953). 
S. Weinberg, Ann.Phys. 194, 336 (1989). 

P. Leifer, On nonlinear nonperturbative modification of quantum mechanics, TAUP Preprint No. 2262 1995(unpublished) 
S.A. Fulling, Aspects of quantum field theory in curved space-time. (Cambridge University Press, Cambridge 1989). 
C. Dewdney and B.J. Hiley, Foundations of Physics 12, 27 (1982). 

C. Dewdney, G. Horton, M.M. Lam, Z. Malik and M. Schmidt, Foundations of Physics 22, 1217 (1992). 



10 



Figure Captions 

1. FIG. 1. Variation of the transmission coefBcient with magnitude V and width a of a scalar potential barrier for 

an incident energy of 1.38. 

2. FIG. 2. Variation of the transmission coefficient with magnitude V and width a of an electrostatic potential 
barrier for an incident energy of 1.38. 

3. FIG. 3. Schrodinger, de Broglie (dotted line) and eigenvector (solid hne) velocities and densities for k = 0.1,a = 
12, V = 0.0306 and R = 0.99, the Schrodinger and de Broglie velocities are not distinguishable, the eigenvector 
velocities are always less than one. 

4. FIG. 4.De BrogHe trajectories and their direction-field for k = 0.1,a = 12, V = 0.0306 and R = 0.99 in the 
region of an interference minima. 

5. FIG. 5. Energy- momentum flow trajectories and their direction- field for k = 0.1,a = 12, V = 0.0306 and 
i? = 0.99 in the region of an interference minima. 

6. FIG. 6.De Broglie (dotted line)and energy-momentum (soHd hne) velocities and densities for k = 0.95, a = 12, V = 
0.36 and R = 0.7. 

7. FIG. 7. De Broglie trajectories and their direction field for k = 0.95,a = 12,V = 0.36 and R = 0.7. 

8. FIG. 8. Eigenvector trajectories and their direction field for k = 0.95, a = 12,V = 0.36 and R = 0.7. 

9. FIG. 9. de Broglie (dotted line)and energy-momentum (solid line) velocities and densities for k = 0.95, a = 
12,V = 4.47 and R = 0.7, in the Klein-region. 

10. FIG. 10. De Broglie trajectories and their direction field for k = 0.95,a = 12,V = 4.47 and R = 0.7, in the Klein 
region. 

11. FIG. 11. Eigen-vector trajectories and their direction field for k = 0.95,o = 12,V = 4.47 and R = 0.7, in the 
Klein region. 
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